Theoretical Evidence for Equivalence between the Ground States of the 
Strong-Coupling BCS Hamiltonian and the Antiferromagnetic Heisenberg Model 
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By explicitly computing wavefunction overlap via exact diagonalization in finite systems, we 
provide evidence indicating that, in the limit of strong coupling, i.e., A/t — » oo, the ground state 
of the Gutzwiller-projected BCS Hamiltonian (accompanied by proper particle-number projection) 
is identical to the exact ground state of the 2D antiferromagnetic Heisenberg model on the square 
lattice. This identity is adiabatically connected to a very high overlap between the ground states of 
the projected BCS Hamiltonian and the t-J model at moderate doping. 
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One of the key questions regarding high Tc supercon- 
ductivity is whether it is possible for a repulsive electron- 
electron interaction alone to give rise to superconduc- 
tivity. Considering that (i) the Hubbard model can be 
transformed into the t-J model in the limit of large on- 
site repulsion, and (ii), at half filling, the t-J model 
reduces to the two-dimensional (2D) S — 1/2 antifer- 
romagnetic Heisenberg model with nearest-neighbor ex- 
change coupling on the square lattice, the answer to the 
above question seems to lie in the relationship between 
2D quantum antifcrromagnetism and superconductivity. 
Although there is a consensus that the ground state of 
the 2D antiferromagnetic Heisenberg model has Neel or- 
der, i.e., staggered spin order yj, the precise nature of the 
ground state wavefunction itself, especially upon doping, 
is controversial. 

The resolution of controversy requires both an unbi- 
ased and unambiguous study. Many analytic approaches 
are based on uncontrolled approximations such as the 
large- N expansion 0. Even numerical approaches such 
as variational Monte Carlo simulation |3, L| are biased 
by the choice of trial state which is assumed to be the 
ground state. Exact diagonalization, on the other hand, 
is a completely unbiased approach. However, its appli- 
cation is limited to finite systems with relatively small 
spatial size. This limitation generates a problem since 
the evidence of an ordered state is usually given by cor- 
relation functions. In the case of superconductivity, the 
relevant correlation function is the pairing correlation 
function: F al3 (r - r') = (cj(r)cj(r + a)c;(r')c T (r' + (3)) 
where a,/3 = x,y. Off-diagonal long-range order can be 
claimed only when F a p remains non-zero in the limit of 
large |r — r'|. Unfortunately, small spatial sizes of finite 
systems accessible via exact diagonalization make the dis- 
tinction between true long-range order and short-range 
order (present even in normal states) ambiguous. 

Despite the small spatial size, however, the size of 
Hilbert space is quite large; depending on the number 
of holes, there are 10 3 -10 5 basis states in the 4x4 lattice 
system (even after translational symmetries are imple- 



mented as reported in this paper) . Since the ground state 
is constructed out of a very large number of possible lin- 
ear combinations of basis states, it is a highly non-trivial 
task to infer a good ansatz wavefunction for the ground 
state itself. Consequently, if there is a method to show di- 
rectly that some form of the BCS wavefunction is identi- 
cal to the exact ground state of the 2D antiferromagnetic 
Heisenberg model, it will provide convincing evidence for 
superconductivity upon doping. In this paper, we put 
forward such a method: we compute the wavefunction 
overlap between the ground states of the t-J model and 
the Gutzwiller-projected BCS Hamiltonian via exact di- 
agonalization. It should be emphasized that, in this pa- 
per, the Gutzwiller projection is applied directly in the 
BCS Hamiltonian instead of being applied onto the BCS 
ground state wavefunction. Note that the Gutzwiller- 
projected BCS wavefunction was originally proposed by 
Anderson Q as a candidate for the ground state either 
on the square lattice with sufficiently strong next-nearest- 
neighbor exchange coupling, or on the triangular lattice. 
The Gutzwiller-projected BCS wavefunction was then 
conjectured to be a good ground state for the t-J model 
on the square lattice at moderate, non-zero doping. The 
difference between the ground state of the Gutzwiller- 
projected BCS Hamiltonian and the Gutzwiller-projected 
BCS ground state is most crucial at half-filling, which will 
be discussed in detail later in this paper. 

The inspiration for using wavefunction overlap comes 
from the fractional quantum Hall effect (FQHE). It is 
well accepted by now that all essential aspects of FQHE 
are explained by the composite fermion (CF) theory [||. 
One of the main reasons why we can put unequivocal 
trust in the CF theory may be the amazing agreement 
between the exact ground state and the CF wavefunction: 
the overlap is practically unity for various short-range in- 
teractions including the Coulomb interaction Q. In this 
paper, we would like to achieve the same methodological 
clarity for the t-J model crucial in establishing the CF 
theory. 

Before we discuss computational details, it is illumi- 
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nating to note that the analogy with FQHE goes much 
deeper than just methodology. To gain a physical in- 
sight into why the CF wavefunction is so accurate, let us 
consider the Laughlin wavefunction |8| which is a sub- 
set of CF wavefunctions at special lowest-Landau-level 
filling factor v = l/(2p + 1) with p an integer; espe- 
cially at v = 1/3, it is given by * 1/3 = J I, ,U, - Zj) 3 
where z = x + iy. Trugman and Kivelson |9j, and 
Haldane showed that ^1/3 is actually the exact 
ground state of a short-range interaction given by V 2 <5(r) . 
For the Coulomb interaction (relevant for experiments), 
the Laughlin state remains extremely close to the exact 
ground state, as shown by practically unity overlap. The 
breakthrough for the CF theory was achieved when Jain 
realized that the Laughlin state is actually composed of 
two parts: the Jastrow factor, ni<j( z i — z j) 2 (dubbed 
flux quanta attachment because of concomitant phase 
winding), and a non- interacting fermionic wavefunction 
of new quasiparicles, i.e., composite fermions. The key 
physical point is that, once the short-range correlation 
is captured by the Jastrow factor, residual correlations 
can be treated as relatively weak long-range correlations 
which are much easier to handle. 

In fact, the separation between short-range and long- 
range correlations can serve as a general method in at- 
tacking strongly-correlated problems. The main question 
is what type of short-range correlation exists in the spe- 
cific problem at hand, and more importantly what func- 
tional form of Jastrow factor describes this short-range 
correlation. Since we are interested in the quantum anti- 
ferromagnetism, it is natural to ask what has been known 
conclusively in the context of short-range correlation in 
quantum antiferromagnetic models. 

To this end, let us consider the ID S = 1/2 quan- 
tum Heisenberg model. This model is important because 
its solution is exactly known for two important cases of 
nearest-neighbor and 1/r 2 exchange coupling. For the 
nearest-neighbor exchange coupling, the Bethe ansatz so- 
lution 11] gives the exact ground state which, for a given 
spin configuration, has an amplitude equal to the product 
of plane- wave states of spin flip excitation with appropri- 
ate phase shift. Since the Bethe ansatz solution is basi- 
cally a product of plane- wave states, it is encouraging to 
guess that some form of Fermi sea state might be closely 
related to the exact Bethe ansatz solution, which turns 
out to be precisely the case. Various numerical works 
0,0] as well as exact analytic studies showed that, 
in addition to closeness in energy, the spin-spin correla- 
tion function of the Gutzwiller-projected Fermi sea state 
has a power-law behavior very similar to the exact re- 
sult. Note that the Gutzwiller projection simply imposes 
the no-double-occupancy constraint. Furthermore, Hal- 
dane and Shastry proved that the Gutzwiller- 
projected Fermi sea state is the exact ground state of 
the ID S = 1/2 Heisenberg model with 1/r 2 exchange 
coupling. 



Now, combined with the fact that the Gutzwiller pro- 
jection is basically an implementation of strong on-site 
repulsion, the above-mentioned similarity leads to a con- 
jecture that the Gutzwiller projection plays the role of a 
Jastrow factor providing the short-range correlation em- 
bedded in quantum antiferromagnetism. To support this, 
we compute the overlap between the ground states of the 
Gutzwiller-projected BCS Hamiltonian and the antiferro- 
magnetic Heisenberg model (in general, the t-J model at 
non-zero doping). In fact, we will show that the ground 
state of the projected BCS Hamiltonian is identical to the 
exact ground state of the Heisenberg model in the limit of 
strong coupling. Also, the overlap is very high (~ 90%) in 
a realistic parameter range relevant to cuprates, which is 
adiabatically connected to the unity overlap in the afore- 
mentioned limit. 

We begin our quantitative analysis by writing the 
Hamiltonian of the t-J model: 

H t -j = V G {H t + Hj)P G , 
H t = -i^(4c Jff +h.c.), 

Hj = jJ2( S i- S J- n i n j/^ (!) 

(i,i> 

where Vg is the Gutzwiller projection operator imposing 
the no-double-occupancy constraint. We obtain the exact 
ground state of H t -j using a modified Lanczos method. 
We have checked that our results for the t-J model com- 
pletely agree with previous numerical studies 0] for all 
available cases. 

Now, let us turn our attention to the Gutzwiller- 
projected BCS Hamiltonian. As mentioned in the be- 
ginning, our approach is rather different from previous 
approaches 0, 0, H 0] which applied the number pro- 
jection as well as the Gutzwiller projection onto the ex- 
plicit BCS wavefunction. We do not take these previ- 
ous approaches for two reasons. First, there is a singu- 
larity due to the spin-density-wave instability inherent 
in the BCS Hamiltonian with on-site repulsion, which 
generates Neel order at half filling. So, it is important 
to study directly the Gutzwiller-projected BCS Hamilto- 
nian instead of applying the Gutzwiller projection onto 
the BCS ground state. The Gutzwiller-projected BCS 
Hamiltonian is given as follows: 

-ffBcs — H t + Ha, 

= A E( c lT c Ui- c a c l*,T+ h - c -) 

i 

~ A E^t c lsa- c !i4 +? M+h.c.), 

i 

^bcs = 'PgHbcs'Pg, (2) 

where H t is given in Eq.Q. Note the sign change in front 
of A in the y direction compared to that of x direction: 
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FIG. 1: Schematic diagram for constructing the ground 
state of the Gutzwiller-projected BCS Hamiltonian, 
4>BCs(Nh,Nh + 2 1 TV), and the exact ground state of 
the t-J model, ipt-j(Nh\N). Nh (N) is the number of holes 
(sites). 



the gap function is given in momentum space by A(k) = 
2A(cos/c;r — cosfcj,). The Gutzwiller projection is built 
in from the onset by working solely in the Hilbert space 
with the no-double-occupancy constraint. 

The second reason is rather subtle, but physically very 
important. Coherent number fluctuations (as opposed 
to incoherent fluctuations in the normal state) are ulti- 
mately responsible for the intrinsic properties of super- 
conductivity. So, coherent number fluctuations should be 
incorporated even into finite-system studies in a funda- 
mental manner. In essence, we diagonalize H^ cs in the 
combined Hilbert space of N e and N e + 2 particles. There 
is, however, a spurious finite-size effect which prevents 
pairing if one naively diagonalizcs H^ cs in the combined 
Hilbert space. In finite systems, the energy cost of adding 
(removing) few particles is not negligible compared to the 
total energy. So, the mixing between states with even a 
few-particle difference is energetically prohibited. We fix 
this problem by adjusting the chemical potential so that 
the kinetic energy plus the chemical potential energy of 
the N e particle ground state is the same as that of the 
N e + 2 particle ground state, which eliminates a spurious 
energy penalty for pairing. Once the chemical potential is 
set, it can be shown that the mixing with other particle- 
number sectors such as the N e + 4 sector is negligibly 
small, even if it is allowed. 

Let us define the following notations: ipBCs(Nh, Nh + 
2 1 TV) denotes the ground state of the Gutzwiller- 
projected BCS Hamiltonian obtained from the combined 
Hilbert space of Nh and Nh + 2 holes in the system of N 
sites. VN h =N denotes the number projection operator 
which projects states onto the Hilbert space of Nq holes 
and re-normalizes the projected states. ipt-j{Nh\N) is 
the exact ground state of the t-J model in the Hilbert 
space of Nh holes in N sites. A schematic diagram is 
shown in Fig 2] 

We now present our exact-diagonalization results of 
overlap for various numbers of holes in the 4x4 square 
lattice system with periodic boundary conditions. Note 
that the 4x4 system is one of the most studied systems 
in numerics because it is accessible via exact diago- 
nalization, yet large enough to contain essential many- 
body correlations. While it is possible to study all pos- 
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FIG. 2: Square of the overlap between the ground state of 
the projected BCS Hamiltonian, VN h =o\tpBCs(Oj 2| 16)}, and 
the exact ground state of the t-J model, |i/) t _j(0|16)), for the 
0-hole case in the 4x4 system (undoped regime). 



sible dopings in the 4x4 system, we concentrate on the 
three most important cases: hole (undoped regime), 2 
holes (optimally doped regime), and 4 holes (overdoped 
regime). Also, we study only an even number of holes 
since, in finite systems, an odd number of holes will ar- 
tificially frustrate pairing order. 

Let us begin with the 0-hole case, i.e., the 2D anti- 
ferromagnetic Heisenberg model. Though the numeri- 
cal evidence for Neel order is quite convincing 0, our 
knowledge of the ground state itself is very limited, con- 
sidering that (i) the semiclassical configuration of stag- 
gered spins is not the exact ground state, and (ii) there 
are still rather strong quantum fluctuations. So, it will 
be satisfactory if one can show that the ground state 
of the Gutzwiller-projected BCS Hamiltonian is a good 
representation of the exact ground state of the 2D an- 
tiferromagnctic Heisenberg model. The criterion for the 
effectiveness of the ground state of the projected BCS 
Hamiltonian is quantified via its overlap with the exact 
ground state of the Heisenberg model. High overlap will 
provide evidence for the existence of pairing which, in 
turn, generates superconductivity upon doping. 

When Anderson [j| proposed the RVB state (which 
is a synonym for the projected BCS state), his insight 
was that electrons are already paired at zero doping, but 
the ground state cannot superconduct (for that matter, 
conduct) because there are no mobile charge carriers. It 
seems reasonable, then, that the ground state becomes 
superconducting as soon as holes are added. However, 
the idea of the RVB state as the ground state of the 
Heisenberg model was rejected because the RVB state 
does not have any long-range ma gnet ic order, while the 
exact ground state has Neel order |2(J, |21| . The situation 
is quite different for the ground state of the Gutzwiller- 
projected BCS Hamiltonian. Fig|2 shows that, at zero 
doping, the ground state of the projected BCS Hamilto- 
nian, 'PAr fe =o|V'Bcs(0j 2|16)}, is actually identical to the 
exact ground state of the t-J model, |-0t-j(O|16)), in the 
limit of A/t — > oo. Since the two states are identical, 
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FIG. 3: Square of the overlap between the ground state of 
the projected BCS Hamiltonian, "Pjv h =2|T/>Bcs(0, 2 1 16) } , and 
the exact ground state of the t-J model, |i/> t _j(2|16)), for the 
2-hole case in the 4x4 system (optimally doped regime). 



it is clear that the ground state of the projected BCS 
Hamiltonian has Neel order. This identity at infinite Alt 
is consistent with previous Monte Carlo simulations |4j, 
in which their variational gap parameter becomes very 
large at small doping. It is important to know that the 
largeness of A/t does not necessarily mean strong super- 
conductivity since, despite strong pairing, there is very 
little charge fluctuation at small doping |22| . 

But, physically, why is the overlap unity at infinite 
A/t, or equivalently t = 0? To answer this, consider 
the t-J Hamiltonian at zero doping, Hj, and the BCS 
Hamiltonian with t = 0, Ha- Hj contains Si ■ Sj which 
prefers the formation of spin singlet pairs. Similarly, Ha 
contains cLcJ, — cj,cj t which also prefers to create spin 
singlet pairs. So, provided that the no-double-occupancy 
constraint is imposed via Gutzwiller projection, Hj and 
Ha should have the same physical effect. 

We now move onto the cases with finite doping. Since 
pairing already exists at zero doping, it is expected that 
the ground state becomes superconducting upon doping. 
We support this by showing that, in the Hilbert space of 
two holes, the ground state of the projected BCS Hamil- 
tonian, VN h =2\ipBcs(^' 2jl6)), has a very high overlap 
(~ 90%) with |-0t_j(2|16)) at optimal A/t for a realistic 
range of J/t: 0.4 < J/t < 0.8 (middle panels in FigEI 
showing the square of the overlap). Note that A/t can 
be taken as a variational paramter. In fact, the optimal 
overlap approaches unity when J/t bocomes sufficiently 
large (bottom, right in Fig[3J). While the large J/t regime 
itself is not very realistic, the high overlap in the realis- 
tic regime is adiabatically connected to the unity overlap 
in the large J/t limit. Incidentally, the identity between 
ipt-j and V'scs at = an d A/t — (top, left in 



Fig|3l is rather trivial because H t .j and H^ cs become 
identical in this case. It is important to note that the 
nature of the identity at large J/t is completely differ- 
ent from that of zero J/t, as manifested by symmetry 
changes of the ground state. The t-J model ground state 
changes its rotational symmetry from s-wave to d-wave at 
J/t ~ 0.08 while the ground state of the projected BCS 
Hamiltonian does so at A/t ~ 0.1. Therefore, the regime 
with large J/t and A/t is completely disconnected from 
the regime with small J/t and A/t. It is important to 
distinguish the rotational symmetry of the ground state 
from the pairing symmetry. The latter is always d-wave 
while the former changes as a function of A/t. 

Finally, we have checked that, in the overdoped regime, 
the ground state of the projected BCS Hamiltonian is no 
longer a good representation of the ground state of the 
t-J model, which is supported by the negligible overlap 
between V Nh =i\^BCS% 4 I 16 )) and IVw(4|16)) for gen- 
eral parameter range. 

In conclusion, we have provided evidence that, in 
the limit of strong coupling, the ground state of the 
Gutzwiller-projected BCS Hamiltonian is equivalent to 
that of the 2D antiferromagnetic Heisenberg model. 
Combined with high overlaps at moderate doping, this 
equivalence supports the existence of superconductivity 
in the t-J model. For future work, it will be interest- 
ing to investigate an analytic approach in proving the 
equivalence. 
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